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ABSTRACT 

We demonstrate that the observationally inferred rapid onset of star formation after parental molec- 
ular clouds have assembled can be achieved by flow-driven cloud formation of atomic gas, using our 
previous three-dimensional numerical simulations. We post-process these simulations to approximate 
CO formation, which allows us to investigate the times at which CO becomes abundant relative to the 
onset of cloud collapse. We find that global gravity in a finite cloud has two crucial effects on cloud 
evolution, (a) Lateral collapse (perpendicular to the flows sweeping up the cloud) leads to rapidly in- 
creasing column densities above the accumulation from the one-dimensional flow. This in turn allows 
fast formation of CO, allowing the molecular cloud to "appear" rapidly, (b) Global gravity is required 
to drive the dense gas to the high pressures necessary to form solar-mass cores, in support of recent 
analytical models of cloud fragmentation. While the clouds still appear "supersonically turbulent" , 
this turbulence is relegated to playing a secondary role, in that it is to some extent a consequence of 
gravitational forces. 

Subject headings: gravitation — instabilities — turbulence — stars:formation — ISMxlouds 



1. INTRODUCTION 

One motivation for pursuing a picture of molecular 
cloud formation by large-scale flows was to resolve the 
"crossing time problem", i.e. the fact that the age 
spreads of young stars are often much shorter than 
lateral dynamical timescales (jBallesteros-Paredes et all 
119991 : iHartmann et all 1200 iL HBBOl). In this picture, 
dense star-forming clouds are produced as flows - ex- 
panding HIT regions, stellar wind bubbles, supernova ex- 
plosions, spiral density waves (e.g. BonncU et al. 200i) 
or gl obal gravitational instabilities (e.g. I Yang et al.l 
- sweep up and condense diffuse interstellar gas. 
The crossing time problem is eliminated because infor- 
mation is not being transmitted along the long dimension 
of the cloud. 

A major attraction of the notion of flow-driven molecu- 
lar cloud formation is that we see several exa mples of it in 
the s olar neighborhood, such as in Cep 0B2 (iPatel et all 
[1998"), in addition to the well-known cloud and star for- 
mation in galactic spiral arms d ue to gas inflow (e.g. 
Elmeereenl 1197 9. 2007; see e.g. iKim etall l2003l and 
Dobbs fc BonnelL 2007b for numerical evidence). While 



it is not always easy to identify specific driving sources in 
all cases, this is not surprising given the complexity ex- 
pected as a result of the interaction of neighboring flows. 
In addition, the presence of molecular clouds well out of 
the galactic plane (e.g., Orion) clearly suggests the need 
for some kind of driving. Given the extensive impact 
that massive stars have on the interstellar medium - HII 
regions, stellar wind impacts, and ultimately supernova 
explosions - it is difficult to see how further creation of 
new star-forming locales by fiows with scales of several 
to tens of pc (or even kpc, in the case of spiral arms) 
could be avoided. 

Star-forming clouds in this picture become somewhat 
accidental associations of gas which are not in virial 
equilibrium, although rough energy equipartition is ex- 
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pected (HBBOl), in adequate agreement with observa- 
tions (Ballesteros-Parcdes 2006). Indeed, it is very diffi- 
cult to prevent global gravitational collapse motions in fi- 
nite c louds of many Jeans masses (jBurkert &: HartmannI 
l2004f ). consistent with the empirical evidence for rapid 
star formation following cloud formation (HBBOl). 

While global gravitational contraction provides 
a plausible mechanism for assembling protocluster 
gas and stars and even overall cloud morphology 
(|Hartmann fc BurkertI l2007t ). it poses difficulties for 
local collapse; non-linear perturbations are probably 
required to avoid swe ep-up in overall collapse modes 
(jBurkert fc HartmannI l2004f ). The necessity of pro- 
ducing star-forming clumps through turb ulent motions 
has long been recognized (jLarso nl Il98l . based on a 
large number o f numerical simulations over the years 
(see review by [Mac Low fc Klcssen 2004). However, 
the source of this turbulence has been a matter of 
debate. The flow-driven cloud formation picture may 
provide an answer, in that the dynamical instabilities 
coupled with rapid cooling and thermal instability which 
naturally result at the shock interface between driving 
material and swept-up gas gener ate turbulence an d 
non-l i near density fluc tuations ( Kovama fc Inutsu 
120021: llnutsuka fc Kovam a 200 | jAudit fc HennebeiL 
12005"; 'Heitsc h et"all l2005t iVazauez-Semadeni et al 
2006; Hcitsch 'et al.l [2b06bl ; fVazquez-Scmaden i et al 
20071: iHenn cbcUc et al. 2007; He nnebelle fc Auditi l20a' 



Heitsch et al.. ,2008b. : Hennebelle et all l2008fl. Based 



on a consideration of timescales. iHeitsch et al.l (|2008aD 
demonstrate that these instabilities do indeed proceed 
much faster than global gravity, as required. 

As simulations of the relevant processes become more 
detailed, we are increasingly in a position to test the 
implications of flow-driven molecular cloud and star for- 
mation. In this p a,per we analyze nume rical models pub- 
fished elsewhere (jHeitsch et al.l I2008U H&08) to illus- 
trate some general predictions of the flow-driven picture, 
and help clarify some confusing issues about molecular 
cloud and star-formation region lifetimes. 
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2. MODELS 

2.1. Overview 

The interstellar medium is filled with flows of vari- 
ous types, many of which result in piling up material 
through shocks. Recognition of this fact is at the heart 
of the flow-driven cloud formation picture. While many 
different types of initial conditions can be envisaged - 
flows driven by hot gas expansion (H II regions, super- 
nova bubbles), flows within mostly molecular gas - we 
focus on the particular case of mildly supersonic flows 
in atomic gas. We do this for two reasons: first, it is 
computationally simplest for our present situation; and 
second, it should be representative of conditions in the 
solar neighborhood, where most of the gas is atomic (or 
ionized) and only a modest fraction in mass and a very 
modest fraction in volume is present as molecular gas. 
We expect that many of the qualitative and even semi- 
quantitative results will apply to other situations, though 
that needs further exploration. 

We also chose to have equal uniform flows entering on 
either sides of our computational volume. This is obvi- 
ously not the most general case. However, we have run 
some models in which differing density atomic flows col- 
lide; this amounts to transforming into the frame of the 
shock front of an expanding bubble (for example). The 
results of these cases do not differ substantially from the 
results we present here, and so we confine the discussion 
to the simplest possible cases. 

We also chose to have smooth inflows and introduce 
only a modest perturbation of their initial interface. We 
did this not out of the absolute conviction that the flows 
themselves are not turbulent, but with the desire of in- 
troducing as little substructure by hand as possible. By 
showing that even relatively smooth conditions rapidly 
introduce turbulence, we can show that any other initial 
substructure will merely add additional structure to our 
clouds. 

Finally, we chose not to use periodic boundary condi- 
tions. This has absolutely crucial effects on the gravi- 
tational accelerations within the forming cloud. Global 
gravitationally-driven flows play important roles in cloud 
evolution, and these cannot be seen in models with peri- 
odic boundary conditions. 

By addressing the formation processes of a cloud, we 
are able to examine the initial conditions which lead to 
structures and velocity fields important for star forma- 
tion. 

2.2. Model Set 

We base our analysis and discussion on the models in- 
troduced in II&08. They describe the flow-driven forma- 
tion of an isolated cloud with the parameters given in 
Table [TJ All models are run on a fixed grid with the in- 
streaming gas fiowing along the x-direction, entering the 
domain (in opposing directions) at the two (y, z)-planes. 
The nominal resolution (i.e. the size of a single grid cell) 
is 8.6 X 10~^ pc. To trigger the fragmentation of the 
(otherwise plane-parallel) region of interaction, we per- 
turb the collision interface. We chose the perturbations 
of the collision interface from a random distribution of 
amplitudes in Fourier space with a top hat distribution 
restricted between wave numbers k = 1..4. 

The inflow is restricted to a cylinder of elliptical cross 



TABLE 1 

Model Parameters 



Name 




Lj^LyLz [pc] 


gravity 


tend [Myr] 


V [pc] 


Hfl 


256 X 5122 


22 X 442 


no 


14.5 


2.2 


Gfl 


256 X 5122 


22 X 442 


yes 


14.5 


2.2 


Gf2 


256 X 5122 


22 X 442 


yes 


14.5 


4.4 



Note. — 1st column: Model name. 2nd column: resolution. 3rd 
column: physical grid size. 4tli column: gravity. 5th column: end 
time of run. 6th column: amplitude of interface displacement. 



section with an ellipticity of 3.3 and a major axis of 80% 
of the (transverse) box size, mimicking two colliding gas 
streams in a more general geometry. IVIodel Hfl is a non- 
gravitating version of Gfl, to compare the role of gravity 
versus that of the thermal instability for the fragmenta- 
tion of the gas streams. 

The inflow density in all models is rij = 3 cm~^ at a 
temperature of Ti = 1800 K and an inflow velocity of 
Vi = 7.9 km s^^, corresponding to a IVIach number of 
M. = 1.5. The flows are initially in thermal equilibrium. 
The models start at time t = with the collision of the 
two flows. The fluid is at rest everywhere except in the 
colliding cylinders. 

To solve the hydrodynamical equations, we used 

the h igher-or der g as-kin et ic grid method P roteus 

JPrendcrgast [l99l ISlvz fc P renderg ai^ [l999|: 

Heitsch et al.ll2004l: iSlvz et al.ll2005l : lHeitsch et al. 2006b!, 
2007f ). The code evolves the Navier-Stokes equations in 
their conservative form to second order in time and space. 
The hydrodynamical quantities are updated in time- 
unsplit form. Self-gravity is implemented as an exter- 
nal source term, also in time-unsplit form. The Poisson- 
equation is solved via a non-periodic Fourier solver, us- 
ing the (IVIPI-parallelized) f f tw (Fastest Fourier Trans- 
form in the West) libraries. The heating and cooling 
rat es are restricted to o ptically thin atomic lines follow- 
ing '^^IfireliOn (|1995'). Numerically, heating and cool- 
ing is implemented iteratively as a source term for the 
internal energy e. While we do not include molecular 
line cooling, the resulting effective equation of state re- 
produces the quasi-isothermal behavior expected at high 
densities, reaching a temperature of « 14 K at densities 
of « 10^ cm~^. 

The x-boundaries are partly defined as inflow- 
boundaries. The inflow is defined within an elliptical 
surface in the (y, z)-plane. The y and z boundaries, - as 
well as the part of the ^-boundaries that is not occupied 
by the inflow -, are open, meaning material is free to 
leave the simulation domain through these boundaries. 
For a detailed discussion of the numerical scheme and 
the models, see H&08. 

In view of the more general situation, we explored a va- 
riety of geometries and flow parameters, such as expand- 
ing shells, unrestricted inflows, and flows with different 
densities, the latter probably being the most general sce- 
nario for e.g. swept-up clouds at the rims of supernova 
bubbles. Despite the variations, the basic mechanisms 
and results do not change substantially. Speciflcally, the 
combination of thermal and dynamical instabilities en- 
hances local gravitational collapse in all cases. 

2.3. CO formation 
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Our models do not explicitly follow the transition from 
atomic to molecular gas. However, conditions in the 
dense post-shock gas should result in molecule forma tion 
at some point (see Glover & Mac Low |2007al l2007bl | for 
a discussion of H2 formation in turbulent flows). To illus- 
trate what would be observed as a molecular (CO) cloud 
we make a simple approximation, motivated by the no- 
tion that CO formation requires shielding by dust grains. 
For each time instance shown in ^j3l we decide whether 
CO is "present" in a particular grid cell by determining 
the expected radiation field integrated over solid angles. 
If the effective UV extinction is equivalent to that of 
an angle-averaged Ay = 1 and the local temperature is 
T < 50 K, we assume CO is present in high abundance. 
Because CO is rapi dly dissociated at lower extinctions 
(jvan Dishoeck fc Black 1988). we do not advect CO for 
simplicity. The radiation field at each grid point is cal- 
culated by measuring the incident radiation for a given 
number of rays and averaging over the resulting sky. The 
ray number is determined such that at a radius corre- 
sponding to half the box size, each resolution element 
is hit by one ray. Thus, fine struc tures and strong den - 
sity variations are resolved (see also lHeitsch et al.ll2006al ) . 
Since our CO maps are post-processed maps, there is 
no timescale of CO formation involved, but CO is as- 
sumed to appear instantaneously in sufficiently shielded 
regions. This simplification leads to "more" CO being 
present, countering the neglection of advected CO. A full 
treatment of CO forma tion including (eve n a simplified) 
chemical network (e.g. iBe rgin et al.ll20d^ . advection of 
chemical species and consistent radiative transfer would 
be a major computational challenge and beyond the de- 
gree of detail necessary for our argument. The issue of 
CO formation is discussed further in § 14.11 

We concentrate on CO rather than on H2 formation 
because CO is the first available tracer to indicate the 
appearance of a molecular cloud, whereas H2 is usually 
no t directly ob s ervabl e in clouds. As Figures 1 and 3 
of iBergin et al.l ()2004D show, H2 formation sets in well 
before CO formation due to self-shielding instead of dust 
shielding - specifically if traces of H2 are already p resent 
in the swept-up material (e.g. Pring le et al.ll2001[ ). 

3. AT THE MOVIES 

We now proceed to look at the evolution of the forming 
cloud with our crude treatment of CO formation. For 
reference, note that in the absence of any substructure 
or collapse, the column density along the x-axis would 
vary with time as 

{Nh) = 2 X 7.9 km X 3 ciR-^t (1) 

= 1.5 X 10^° cm"2(i/Myr) , 

and the mass of the cloud would be 

Me ~ 5OOM0(t/Myr) , (2) 

using the approximate area of our model cloud. Thus, 
at the end of the three simulations (t = 14.5 Myr), the 
average column density through the cloud in the absence 
of substructure would be Nh ^ 2.2 x 10^^ cm~^ or Ay ~ 
1.1, and the cloud mass Mc ~ 7500Mq. The actual 
molecular mass of the cloud will be smaller than this 
(see below), so our simulation is related to the formation 
of a small molecular cloud. 



3.1. Cloud Formation and Onset of Star Formation 

Figures [T] and [5] provide snapshots of the cloud evolu- 
tion at various times in the form of column density maps, 
where we use contours to illustrate H I levels at levels of 
log Nhi — 20, 21, 22, with colors showing those hydrogen 
column densities at which CO is "strongly present" . Fig- 
ure [1] shows the cloud formation seen along the inflows, 
and Figure [2] gives an "edge-on" view. 

While at early times the small-scale structure in the 
HI column densities indicates the rapid fragmentation of 
the inflows due to thermal and dynamical instabilities in 
all three cases, the subsequent evolution is very different. 
Consider first the left-most columns in Figure [1] which 
illustrates the evolution of model Hfl, in which grav- 
ity has been turned off, seen along the direction of the 
flows. Very little CO formation occurs over the ^ 15 Myr 
buildup of material. The resulting "cloudlets" are small, 
relatively thin (Ay ~ 1), and in some cases transient, 
due to the continuing input of turbulent energy from the 
inflows. 

Things are dramatically different for the evolutionary 
sequence Gfl in the middle column of Figure [H where 
gravity is present at all times. Although there is lit- 
tle change in behavior a,t t = 5.3 Myr, by t = 7.6 Myr 
Gfl has significantly more "CO" , although the concen- 
trations are not strongly self-gravitating f H4.2p . By 
t = 9.9 Myr the model with gravity has formed a bro- 
ken ring. This structure results from the "edge effect" 
of non-linear gravitational a cceleration in a finite cloud 
(|Burkert fc HartmannI [200I HfcOS); it IS m some sense 
an "echo" of the initial cloud boundary. The slightly 
more compressed appearance of the HI contours (see 
also Fig. 2 of HfcOS) indicate that it is the global grav- 
itational modes and the resulting lateral collapse mo- 
tions (perpendicular to the inflow), which increase the 
local column density. The peak column densities at 
10 Myr are above 10^^ cm~^, compared to the aver- 
age value ^ 1.5 x 10^^ cm~^ expected at that time. 
Overall, the rapid production of shielded molecular gas 
is strongly enhanced by lateral gravitational collapse or 
sweep- up of gas, an effect initially suggested by Bergin 
et al. (|2004l : see UM- The lateral collapse overcomes 
the perceived limitations of the flow-driven cloud for- 
mation scenario that within reasonable timescales the 
column densitie s would be too low for ef ficient molecule 
formation (e.g. iMcKee fc Ostrikeill2007l ). This can also 
be seen in Figure [31 which shows the density-weighted 
width of clouds Hfl and Gfl along their short axis against 
time. At around 8 Myr, the widths start to separate, 
indicating the global gravity effects. The effect in the 
widths is somewhat lessened by the edge effect due to the 
non-linear gravitational accelerations, which nevertheless 
help sweeping up material laterally and thus increasing 
the shielding. By t = 12.2 Myr (Fig.[T]) several local small 
condensations appear which are very dense and strongly- 
self gravitating. The overall region continues to collapse 
laterally under gravity until it forms something close to 
a single filament at t = 14.5 Myr. 

The evolution of model Gf2 in the right-hand column is 
fairly similar to that of Gfl. The larger physical pertur- 
bations result in somewhat denser concentrations of gas 
by t = 9.9 Myr, and the overall collapse tends to produce 
a single narrow filament by the end of the simulation. 
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Figure [2] shows column densities in a side view of the 
three simulations. The main concentrations formed in 
models Gfl and Gf2 are driven initially at the positions of 
inflection s in the initia.1 pert urbed interface by dynamical 
focusing (Hucckstacdt 2003; H&08) due to the non-linear 
thin shell instability ( Vishniac.1994) . These views show 
that the resulting clouds and dense concentrations are 
far from spherical or even spheroidal, but quite elongated 
and, eventually, strongly filamentary. 

Once column densities oi N 10^^ cm~^ have been 
reached (at « 10 Myr), gravitationally dominated cores 
start to form (symbols in Fig. |4]). Core s are identified 
by a modified versio n of CLUMPFIND (jWilliams et all 
119941: iKlessen et all |2000[ ). with the additional condi- 
tion that the thermal plus kinetic energy content is 
less than half of the gravitational potential energy (see 
H fc08). This is consistent, for example, with the results 
of lOnishi et"aI1 (|1998D . who studied cores in Taurus with 
C^^O and found a column density threshold for star for- 
mation of about 8 X 10^^ cm~^. 

We cannot follow the collapse of these dense regions 
to anything approaching the sizes of stars. However, it 
is instructive to consider not only whether the regions 
are self-gravitating but the timescales upon which they 
might collapse. Figure [5] shows the distribution of mass 
as a function of its local free-fall timescale, binned lin- 
early for the top row and logarithmically for the bottom 
row for clarity. The bulk of the mass in the simulation 
domain remains at long free-fall times. The formation of 
the first gravitationally-bound cores occurs in the simu- 
lations with gravity at about 10 Myr, as shown in Figure 
131 at about the same time, a small amount of mass ex- 
hibits free-fall timescales of order 1 Myr or so. It is only 
ati'^ll — 12 Myr that appreciable amounts of mass 
are present with free-fall timescales of < 1 Myr; thus, we 
estimate that this is the epoch at which star formation 
begins. It is evident, given the increase in the number of 
self-gravitating cores and the decrease in free-fall times, 
that star formation in this model would be an increasing 
function of time to the end of the simulation. 

4. DISCUSSION 

4.1. What would we observe? 

The cloud evolution that we have presented obviously 
takes far too long to be observable. For purposes of com- 
parison with observation one can think in terms of ob- 
serving four different clouds, each having the same evo- 
lution, but studied at different epochs. We also compare 
at each stage with some observed clouds of similar prop- 
erties, with the primary goal of checking whether our 
prescription for CO formation appears reasonable, either 
with respect to observations or to theoretical treatments 
for these clouds. We restrict attention to model Gfl ini- 
tially. 

"Cloud 1" (the top row of Figure [1]) would be difficult 
to detect in molecular gas emission, but could be inferred 
through absorption line studies and extinction, especially 
in the small blobs. This object would be classified as a 
diffuse H I cloud with a visual extinction Ay ~ 0.4. 

"Cloud 2" (second row), with an age of 7.6 Myr, would 
still be classified as a diffuse H I cloud. It has a small 
amount of CO, mostly in localized positions; Figure |4] 
demonstrates that very little of the cloud is CO-bearing 
at this phase. The average extinction through the cloud 



would be Ay ~ 0.6, with of course some localized higher 
column density regions; such properties are roughly con- 
sistent with local lines of sight to, for example, C Oph and 
C Per, for which van Dishoeck & Black's (1988) models 
are consistent with the presence of small amounts of CO. 

Here we should emphasize that our prescription for CO 
formation relies heavily on results such as those of Fig- 
ure 5a of Ivan Dishoeck fc Blackl (|1988f ). which indicate 
photodissociation timescales > 1 Myr for Ay > 1. Their 
analysis implies the presence of substantial amounts of 
H2 - roughly comparable to the amount of H in the 
( Per and ^ Oph lines of sight. However, our models 
would not form any significant amounts of H2 at the 
stage of "Cloud 2", because essentially none of the gas 
is at densities hig h enough t o form molecular hydrogen 
rapidly enough (Bergin et al.i 2004; formation timescales 
> 1 Myr at densities < 10'^ cm~'^, corresponding to local 
free-fall timescales > 1 Myr; see Figure [5]). The only way 
that there would be significant amounts of H2 present at 
this time (or in " Cloud 2") is if it we re present in the in- 
flowing gas (e.g.. iPringle et ani2001[ ). Thus, if anything 
we may have overestimated the amount of CO present at 
this stage. 

"Cloud 3" (third row) now begins to show significant 
amounts of CO. According to FigurelH the mass in CO is 
now about SOOA/q in the small-perturbation model Gfl 
and about 10"^ M© in the large-perturbation model Gf2. 
However, no self-gravitating cores have formed at this 
point. 

One might regard "Cloud 3" as an example of 
the most massive of the so-called translucent clouds, 
which are generally observed a t high galactic lati- 
tude (e.g., iMagnani et al.l 119851 : iMagnani et~an Il996t 
lYamamoto et al.ll2003| ). The constraint on galactic lati- 
tude probably represents an observational selection effect 
rather than a complete absence of such objects in the 
galactic plane. Observationally, star formation is rare or 
absent in this type of cloud, consistent with our simu- 
lations. For example, in their study of the HI filament 
associated with the molecular complex including MBM 
53, 54, and 55, which has a mass of about I2OOM0 and 
a length of about 30 pc (assuming a distance of 150 pc), 
lYamamoto et al.l ()2003l ) found no evidence for any star- 
forming cores, though there is evidence for one or two 
young weak-emission T Tauri stars in the region. Note 
also that translucent clouds are known to have variations 
in CO abundances of an order of magn itude, even for 
vario us positions within the same cloud ^Magnani et al.l 
Il998l and references therein) which of course is what 
would be expected for a turbulent, structured translu- 
cent cloud like "Cloud 3" . 

"Cloud 4", which appears at i 12 Myr, is now a 
molecular cloud, with a mass associated with CO of order 
1200 Mq. It now contains a small but significant amount 
of mass which is both gravitationally bound and is at 
densities corresponding to free-fall timescales less than 
1 Myr; thus we consider that it has now begun to form 
(low-mass) stars. 

Finally, "Cloud 5" is now a full-fledged molecular 
cloud, with a mass associated with CO of order 2500 
Mq and dense cores comprising of order at most 25% of 
the mass associated with CO and of order 5% of the total 
mass of the cloud including H I. The cores are contained 
within a strongly confined region spatially, due to the 
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Fig. 1. — Time sequence of column density maps for models Hfl (left), Gfl (center) and Gf2 (right) seen along the inflow. Each frame 
measures 22 X 22 pc^. Contours show HI column density, at log Nhi = 20, 21, 22, and colors show column density of gas containing CO in 
the range as indicated. Self-gravity enhances the amount of CO. 



6 



Rapid Cloud and Star Formation 



tm.a 



I21.S 




Hfl, t= 5.3Uff 



I3J.5 



120.0 



1 ^" 



I S3 
21.5 
30.0 




.J- 



Hfi, t= s.'si^ji:^ » 



>■ 



Hfi. l< 



ISl.5 




lai.5 



Gfi, t= 5,awyij 



\ 1' 



lEO.O 



131.5 





j£O.C 



133.0 '.IvSi. . 



121.5 



.120.0 



\ 4 
(1, 1=14,5^ 




IS3.Q 



121.5 



!20,0 



133,0 JV?^tf-< 



.1 eo.o- 




Fig. 2. — Time sequence of column density maps for models Hfi (left), Gfl (center) and Gf2 (right), seen perpendicular to the inflow. 
Each frame measures 11 X 22 pc^. Contours show HI column density, at logNui = 20,21,22, and colors show column density of gas 
containing CO in the range as indicated. 
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Fig. 3. — Mass profile width along the shorter (vertical in Fig.[TJ 
axis of model clouds Hfl (solid line) and Gfl (dashed line) against 
time. Global collapse leads to a narrower width for model Gfl. 
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Fig. 4. — Cloud masses in H I (thin lines), CO (thick lines), and 
in all gravitationally bound cores (symbols) against time for models 
Hfl, Gfl and Gf2. Approximately 2 Myr after CO formation, the 
first cores form (see text for discussion). The turnover in the H I 
and CO masses for model Gfl is an artifact caused by switching 
off the cooling for densities n > 10^ cm""^ (see Fig. 11 of H&08 for 
a detailed explanation). 

global gravitational collapse of the material. 

The behavior of model Gf2, with a larger initial per- 
turbation, is very similar, except that high densities are 
achieved a little earlier, the cloud collapses into a better- 
defined filament, and the efficiency of massive core/star 
formation is larger (Figured]). 

Note that the fraction of mass in gravitationally-bound 
regions does not directly correspond to the rate of star 
formation. For example, in model Gf2 at the end of the 
simulation there is about 2OOOM0 in gas containing CO, 
but only ~ IO'^Mq of this material is actually contained 
in regions with free-fall times of 1 Myr or less (Figure [5|). 
Thus in all cases we expect star formation efficiencies 
would be on the order of a few percent at most, by the 
end of the simulations. 

The sizes and masses of models Gfl and Gf2 at the 
final, "Cloud 4" stage, are reasonably consistent with 
those of on e of the major fil aments in the Taurus molec- 
ular cloud ([Hartmannl 12002) . One may even take this a 
step further and argue that star (core) formation could 
have only started in these models at times ^ 11 — 12 Myr, 



so that at 14.5 Myr the age spread is only of order 2-3 
Myr, again reaso nably consistent with that of Taurus 
(jHartmannl |2003[ ) : but one would then need to shut off 
star formation to maintain the agreement with observa- 
tion (§4.6). 

The above sequence of "clouds" is not meant to suggest 
that all low-mass atomic clouds will eventually become 
molecular, star-forming regions. In particular, high- 
latitude clouds may well remain atomic, simply because 
they cannot sweep up enough mass. 

4.2. Rapid molecular cloud formation 

The simulation sequences shown in Figures [T] and [2] 
support previous suggestions of rapid cloud and star for- 
mation. Specifically, HBBOl argued that the paucity 
of substantial molecular clouds without star formation 
meant that collapse of molecular cloud cores must fol- 
low soon after cloud formation. They further suggested 
that this was the result of the required shielding column 
density for molecular gas against the dissociating inter- 
stellar radiation field, which was similar to that needed 
for gravitational collapse. Our simulations demonstrate 
this behavior, with molecular (CO) gas setting in near 
t ^ 10 My and core formation soon after, with filamen- 
tary gas structure. 

A few principles explain why our simulations show the 
observationally required behavior. One is that the cloud 
is formed with relatively high densities even in the atomic 
phase, and thus is "pre-formed" with a short free-fall 
time making rapid collapse possible. Figure [6] shows that 
the pressure profiles along the infiow axis, averaged over 
the surface area of the cloud for model Gfl, roughly re- 
flect the ram pressures of the inflows until gravity be- 
comes important. This means that the density within 
the cloud will approach 

p Piv1{kBT/iimH)~^ , (3) 

where pi and Vi are the inflow density and velocity, re- 
spectively, and T is the temperature. In terms of densi- 
ties and velocities pi and Vi relative to the values used in 
these simulations (u,; ~ 7.9 km s~^ and = 3 cm~^, see 

p ^ 1.3 X IO-21T3-1 p.f}fgcm-3 , (4) 

or 

n - 750T3o^ p,if cm-3 , (5) 

where Tso is the temperature measured in units of 30 K, 
typical of the cold neutral material. (Note that in the 
simulations the mean molecular weight is assumed to be 
unity, as used in these equations.) The free-fall time is 
then 

-"^(3i;)"'-''»^^^'-'"*.-'My.-. (6) 

This estimate is in reasonable agreement with Figure [SJ 
which shows a substantial amount of gas at r^y ^2 Myr, 
with higher temperature material at r/j ~ 6 Myr. Thus, 
simply due to ram pressure, moderately cold gas in the 
cloud is formed at densities such that free-fall timescales 
are short. This allows rapid gravitational collapse once 
sufficient mass has been accumulated. 

Ano ther feature of our model is that, as iBergin et al.l 
(|2004f) showed using a one-dimensional shock model with 
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Fig. 5. — Time evolution of mass at given freefall time, for models Hfl, Gfl and Gf2. Top: Linearly binned free-fall time, bottom: 
logarithmically binned freefall time. The bulk of the mass stays at freefall times longer than the simulation time. At t s:^ 10 Myr, a few 
100 Mq have reached freefall times Tff < 1 Myr, indicating local collapse. Even when logarithmically binned, the mass at small freefall 
times is sufficient to form stars. 



chemistry, much of the time spent accumulating cloud 
mass is spent in the atomic phase; such clouds would not 
be recognized as "proto" molecular clouds. Beyond this, 
however, local gravitational collapse tends to occur with 
CO formation because global gravitational collapse is re- 
sponsible for accelerating the formation of high-column 
density, highly-shielded gas which can become molecular. 
This was already foreseen by Bergin et al., who suggested 
that lateral collapse under gravity would be responsible 
for producing "runaway" increases in column density, en- 
hancing the rate of molecule formation beyond what can 
be achieved in a one-dimensional flow. We observe pre- 
cisely that behavior in our cloud sequence in Figures 
[1] and [2] (compare the left-hand column without grav- 
ity to the ce nter and right-hand columns with gravity). 
iDobbs et all ^08) discuss the formation of H2 clouds in 
spiral arms and argue that substantial fractions of H2 can 
be formed rapidly even without self-gravity. This does 
not contradict our findings, since - as discussed in H2.3I - 
H2 is expected to form well before CO, at lower column 
densities. In addition, the densities of the clouds foimd in 
these large-scale simulations are generally smaller than 
typical molecular clouds in the solar neighborhood. Fi- 
nally, even if H2 can form without self-gravity, this does 
not mean that both H2 and CO will not form faster with 

s elf-gravity. 

iBurkert fc Hart mannI ()2004f ) showed that cold sheets 
were also susceptible to global gravitational collapse. Ig- 
noring concentrations due to edge effects, which can oc- 
cur much faster, they found that the timescale for global 
sheet collapse with outer radius r is 



1/2 



(7) 



If we ignore the collapse motions and simply allow mass 
addition at a rate such that E — piVit, then we can 
estimate the distance scales over which global collapse 



should be present as 

r - 5.4t^p,{i,pc, (8) 
where ^7 is the time measured in units of 10 Myr. Thus 
it is not surprising that global gravitational collapse is 
strongly underway at the time of massive core condensa- 
ti on in the simulations. 

IBurkert fc HartmannI ()2004l ) also demonstrated that fi- 
nite sheets experience highly non-linear gravitational ac- 
celeration as a function of position which tends to cause 
material to pile up near cloud edges. Our simulations 
confirm this; even though the resulting clouds are three 
dimensional instead of being sheets, they are flat enough 
that the initial elliptical boundary is echoed in the result- 
ing fllamentary collapse (an e ven stronger versio n of the 
edge effect is present in lVazauez-Semadeni et al.l (j2007l) . 
whereas we took steps to suppress the strength of the 
edge effect by reducing the inflow mass addition near the 
boundaries). 

To summarize, reasonable ram pressures of inflowing 
material result in cold atomic gas with densities suf- 
ficiently large that gravitational collapse can be rapid 
once sufficient material is accumulated. Substantial tur- 
bulent fragmentation is needed to form stellar mass ob- 
jects. Global collapse is generally important in driving 
up column densities and enhancing the formation of fila- 
mentary structures through edge effects. Free-fall times, 
Jeans masses and lengths will decrease, and global col- 
lapse will occur faster for higher inflow rates and ram 
pressures, which probably are needed to make more mas- 
sive clouds than those of our simulations. 

We emphasize that simulations with periodic bound- 
ary conditions cannot capture the evolution seen here, 
in particular the rapid increase in column density within 
filamentary structure. Simulations with periodic gravity 
can only be relevant if started with substantial perturba- 
tions, and limited to scales much smaller than the overall 
sizes of clouds. 
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Fig. 6. — Pressure profiles along the inflow axis for model Gfl. 
Shown are total, dynamical ("kin"), thermal ("int") and gravita- 
tional ("grv") pressure. 



4.3. "Accelerated" Star Formation 

Palla & Stabler (|2000l . |2002[ ) estimated ages for stars 
in several local star-forming regions, finding that the rate 
of star formation generally showed an increase up to the 
last 1 Myr or so. As Hartmann (|2002l l2003f ) pointed 
out, it is implausible that all of these star- forming re- 
gions should be coordinated in their evolution; it must 
be that the typical lifetimes of star-form ing regions are at 
most a few Myr. Hartmann (|2002L[2003i ) also pointed out 
some problems related to neglect of observational errors, 
logarithmic vs. linear binning in age, and likely isochrone 
or birthline problems which exaggerate the age spreads 
of these regions. Nevertheless, there is evidence that at 
a small fraction of the stars associated with these areas 
are a few Myr older than the bulk of the population. 

Accelerated star formation is expected in a trivial 
sense, since starting from a rate of zero star formation 
to a finite star formation rate mathematically requires 
an acceleration. However, it also occurs naturally in our 
evolutionary model of molecular cloud formation. The 
Gfl and Gf2 simulations shown in Figures [T] and [5] show 
that a few dense concentrations show up by 10 Myr or so 
(see also the core mass history. Fig. 5 in H&08). One can 
easily infer that, if we had sufficient resolution, we would 
observe a few stars being formed from the very densest 



of the dense concentrations first; later on, as global and 
local collapse proceeds, more and more stars would form. 

Our simulations emphasize that the global age spread 
in a star-forming region is an upper bound to the 
timescales of local collapse, as compressions will never 
be perfectly coordinated in time across any given cloud. 

4.4. The Onset of Star Formation 

In our model, turbulence generated through the cloud 
formation process provides a mechanism to produce ther- 
mal and dynamical fragmentation. While this fragmen- 
tation is highly efficient in generating cold high-density 
cloudlets, the "pre-formation" of such cloudlets is lim- 
ited by the minimum Jeans mass achievable. With Cs 
denoting the (isothermal) sound speed, we have 



M.7 = 



G 



3/2 



,-1/2 



(9) 



where Pram — Pi'i'i i again in units of our simulations. To 
form e.g. solar- mass objects, higher pressures are needed. 
Since invoking "supersonic turbulence" as a fragmenta- 
tion mechanism in its own right without specifying its 
physical source is somewhat unsatisfactory, the only re- 
maining free energy source leading to a self-consistent 
pict ure of the f ormation of prc-stcUar cores is gravity 
(e.g. iBertoldi fc McKce 1992; Field ct al. 2008). 

We have already seen that gravity helps along the "CO 
formation" in our clouds (Figs. [1] and [2]), and also that 
unless the simulations include gravity, they are not go- 
ing to form the high-density cores with freefall times 
Tff < 1 Myr (Fig. [5|). Figure [7| demonstrates that it 
is gravity which eventually leads to the high pressures 
required for fragmentation into small masses. It shows 
the mass-weighted pressure histograms of models Hfl 
and Gfl in the same time sequence as in e.g. Figure [TJ 
At t — 5.3 Myr, the distributions arc indistinguishable. 
With increasing time, gravity drives more and more mass 
to higher pressures, until the pressure distribution follows 
the Jeans mass (eq. [5]) at the lowest temperature thresh- 
old (indicated by the blue dashed line). In contrast to 
model Gfl, the pressure distribution of model Hfl stays 
more or less constant with time. 

4.5. Turbulent Support 

It has been claimed tha t the loss of turbulent support 
(e.g. iHuff fc Stahleil[2007b would l ead to an (accelerated) 
star f ormation in molecular clouds (jPalla fc Stahlei|[2000l 
I2002f) . A closer look at the energy budget (Fig. [S]) in our 
model clouds however reveals that in the context of fiow- 
driven cloud formation, it is actually the build-up of mass 
and thus the deepening of the gravitational potential well 
which leads t o collapse, but not the loss of tu rbulence (see 
also Fig. 8 of lVazquez-Semadeni et al]l2007f ) . The energy 
densities were calculated using only gas at T < 50 K, ap- 
proximately corresponding to the molecular (Gfl, Gf2), 
but at least cold, dense gas (Hfl). The gravitational ener- 
gies grow faster with time than the kinetic energies drop - 
and eventually, once the molecular cloud has "appeared" 
(see Figs. [Hand [5]), the kinetic energies start to rise, in- 
dicating that the turbulent "motions" at that stage are 
driven to some extent by gravitational collapse (see also 
lBallesteros-Paredesl2006t IVazguez-Semadeni et al.ll2007l : 
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Fig. 7. — Mass (in Mq) within logarithmic pressure bins for 
models Hfl and Gfl, at the same times as in e.g. Figure [l] Colors 
indicate temperature thresholds. The left vertical black dashed line 
stands for the thermal pressure of the inflows, while its right-hand 
counterpart shows the ram pressure of the inflows. The slanted 
lines in the temperature color code indicate the Jeans mass (eq. [9]l. 
With increasing time, gravity drives more and more mass to higher 
pressures, allo wing the fragmentation into successively s maller 
structures ifFleld et al.ll200g : TVazquez-Semadeni et aLllSCM . 

iField et al.ll200a IVazguez-Semadeni et ai]|200 8' for a re- 
cent discussion). In other words, while the turbulence 
might decay slightly during the formation of the model 
cloud (in the atomic phase), it is increasing - driven by 
gravity - during the molecular phase. One would ex- 
pect stellar feedback to even strengthen this trend, if 
anything. Figure [5] demonstrates that despite the gravi- 
tational dominance, the linewidths still can look "turbu- 
lent". The linewidths shown were calculated for all gas 
at T < 50 K, along a single line-of-sight of 16 resolution 
elements width through the center of the cloud. 

We note that the kinetic energies shown in Figure [8] 
are subvirial by a factor of up to 3 at late times, and 
that the linewidths (Fig. ^ are on the smallish side. We 
believe this is due to the combination of three issues. 
First, the head-on colliding flows do not allow for large- 
scale shearing motions. Those would be expected e.g. 
in molecular clouds formin g in spiral arms . Numerical 
models of this process (Bonnell et al.lf2b06f ) including a 
shear c omponent rep roduce turbulent linewidth-size re- 
lations (jLarso I3 fT98lh . Second, non-uniform inflows will 
increase the level of turbulence in the resul ting cloud, as 
demonstrated by Dobbs & Bonncll ( 2007a[ ). They found 
the slope of the lincwidth-sizc relation to flatten with in- 
creasing filling factors of the clumpy inflows. Third, the 
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Fig. 8. — Energy densities (or pressures) of the three isolated 
cloud models Hfl, Gfl and Gf2, for gas at T < 50 K, approximately 
corresponding to the molecular gas. The internal energy (black 
lines) stays constant, while the gravitational energy (green lines) 
increases with time, offsetting the slight decrease in kinetic energy 
(red lines). 

dynamical range for gravitational collapse in our mod- 
els is somewhat limited, possibly suppressing turbulent 
motions once the cloud has collapsed. 

We take the discussion of the perceived loss of turbu- 
lent support as an opportunity to demonstrate how the 
star formation in our models (strictly speaking, we are 
only forming cores) would b e represented by the normal - 
ized "star forma tion rat e" of [Krumholz fc McKed ()2005D 
and lKrumholTfc Tan ( 20071). 

M 

SFRff = — -Tff. (10) 

M(>n) ^ ' 

Here, M(> n) is the mass of the cloud above a density 
threshold n, M is the star formation rate, and r// is 
the freefall time at a given density. The choice of this 
density obviously will affect the results. Figure [TOl shows 
the SFRff as defined above, for various times against 
the threshold density n (model Gfl). As a proxy for the 
star formation rate M we use the time derivative of the 
mass accretion history (see e.g. symbols in Fig. 2]). This 
includes the formation of new cores as well as the mass 
accretion onto already formed cores. It is also consistent 
wi th the est imates of M fr om the simulations discussed 
bv lKrumholz fc Ta^ (|200l . 

At early times {t = 11.4,12.2 Myr), the models are 
consistent with the ana lytical predictions (see Fig. 5 of 
iKrumholz fc TanI [2007f) : the SFRff ranges around a 
few percent. Obviously, it is a matter of interpretation 
whether this low percentage is to be seen as "slow" star 
formation: especially at the lower density thresholds, the 
freefall time Tff is not representative of the local con- 
ditions under which the "stars" form. At later times, 
several effects are discernible. For low density thresh- 
olds (i.e. considering the "whole" cloud, see also Fig. [5]), 
still SFRff < 10%. For higher density thresholds (i.e. 
selecting for the regions that are locally collapsing in 
the simulation) , the SFR f f increases but stays short of 
« 50%. This is not surprising, since (a) these regions are 
not resolved, i.e. they would not convert all their mass 
into stars so that the "real" SFRff would be substan- 
tially lower and (b) our models do not include feedback. 
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Fig. 9. — Line-of-sight velocity distribution for models HFl (top) and Gfl (bottom), for a single line of sight parallel to the shorter (left) 
and longer (right) axis of the cloud perpendicular to the inflow. Three times as in Figure [T] 
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Fig. 10. — The sta r formation rate normalized by the free fall 
time SFRff (eq. 1101 ). The plot has been generated following the 
recipe given bv IKrumholz fc TanI I I2007I ). i.e. the free-fall time is 
computed for the mean (left) and median (right) density above 
the density threshold given on the a;-axis. Symbol sizes indicate 
simulation times. 

Moreover, the models discussed here are an extreme real- 
ization of a range of scenarios whose other extreme would 
be a perfect shear flow (Williams et al., in preparation). 
With a larger component of shear or angular momentum, 
the SFRff could be reduced even further. 

4.6. Other limitations 

In addition to limited resolution, our models do not in- 
clude stellar energy input, which is the the most plausible 



mechanism to reverse the collapse of our self-gravitating 
molecular clouds. It is difficult to see how stellar en- 
ergy input can exactly balance the non-linear variation 
of gravity with position within the cloud; there is more 
parameter space for either expansion or contraction. We 
note that even at the end of our simulations, most of 
the mass of the cloud, including that of the molecu- 
lar (CO) regions, has long free-fall times and low den- 
sities. The low filling factors indicated by the SFRff at 
n « 100 cm~^ and by Figure [5] suggest that stellar feed- 
back can in principle blow away most of the cloud mass 
and thus keep overall star formation efficiencies low (see 
also Elmcgrcen 2000) . While this remains an important 
issue of our scenario which needs further investigation, it 
should also be pointed out that to our knowledge there 
are no simulations of long-lived clouds that are not com- 
puted with periodic boundary conditions, the adoption 
of which begs the question. 

Magnetic fields coul d in princi ple suppress or affect 
the fr agmentation (e.g. | Fieldlll965i for the thermal insta- 
bility; iHeitsch et al.l 120071 for the dynamical instabilities 
in colliding flows), although this capability seems to de- 
pend strongly on the number of dimensions considered 
(e.g. Ilnoue fc Inutsu ka 2008). Also, the sheer accumula- 
tion of mass along fieldlines might render the fields less 
important than com monly thought, as demo nstrated in 
numerical models bv iHennebelle et al.l (|2008f ). 
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5. SUMMARY 

Making a simple but reasonable approximation for the 
formation of CO from atomic gas, we have shown that the 
formation of molecular clouds from large-scale flows can 
be tied closely to the formation of self-gravitating struc- 
tures. Our numerical simulations support many parts 
of the rapid cloud and star formation scenario, in par- 
ticular the rapid appearance of CO in self-gravitating 
clouds, quickly followed by the formation of dense cores 
with free-fall timescales of order 1 Myr, and that global 
gravitational collapse plays an important role in cloud 
evolution. The simulations reproduce expected low star 
formation rates at early times. For the long-term evo- 
lution of the star formation history and the molecular 
cloud itself, the issue of stellar energy input (feedback) 



needs to be addressed to develop a more comprehensive 
picture of star formation in the solar neighborhood. 
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